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1.  INTRODUCTION 

Thennally-induced  error,  as  the  term  is  qjplied  to  machine  tools,  is  that  error 
induced  in  the  cutthtg  or  forming  process  as  a  r^ult  of  temperature  changes  to  and  the 
resulting  growth  in  the  structure  of  the  machine.  The  heat  disnpated  from  motors  and 
beatings  and  other  heat  generating  components  is  almost  always  unevenly  distributed 
within  the  machine  structure.  This  causes  structural  deformation  that  is  referred  to  as 
thermal  growth.  Thermal  growth  is  often  the  miyor  cause  of  machining  errors  in  computer 
numerical^  controlled  (CNC)  machining  centers.  In  some  cases,  thermal  growth  can 
contribute  up  to  70%  of  the  total  machining  error. 

As  the  demand  has  increased  in  recent  years  for  better  quality  control  and  tighter 
manufimturing  tolerances,  more  research  efforts  have  been  directed  to  the  understandmg 
and  control  of  the  thennal  growth  of  CNC  machining  centers.  Although  various 
techniques  fr>r  monitoring,  characterizing  and  controlling  spindle  thermal  growth  have 
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been  devdoped  and  d<»nonstrated  in  a  number  of  universitiM  and  research  institutes,  the 
bulk  of  the  work  is  still  too  costly  and  cuiiU)ersome  to  be  practical  in  a  production 
environment. 

Several  problems  are  common  for  engineers  and  researchers  working  with  thermal 
growth  control:  (1)  a  lack  of  standard  equipment  to  perform  quality  measuremoits  on 
machines,  (2)  a  lade  of  standard  user-fiiendy  software  and  amdysis  tools  for  machine 
characterization,  (3)  too  much  time  required  in  the  characterization  and  nuxleling  process, 
and  (4)  the  model  is  generally  inadequate  for  dealing  with  the  situations  encountered  in 
actual  production  processes  and  environments. 

In  June  1991,  Automated  Predsion,  Inc.  (API)  obtained  a  research  and 
development  contract  from  the  U.S.  Air  Force,  Wright  Laboratory,  Manufttctuiing 
Technology  IMrectorate  to  develop  a  commercially  viable,  robust  and  user-friendlv 
spindle  thermal  growth  characterization  and  compensation  system  for  production  use  on 
CNC  machining  centers.  Acting  as  a  collaborator  and  potential  end-user  of  the  system. 
Monarch  Machine  Company  in  Cortland,  New  York  joined  the  program  by  providing  test 
&dliti^  CNC  machining  centers  and  en^neering  time.  The  entire  program  consisted  of 
the  following  major  tasks: 

(1)  development  of  a  low-cost,  high-performance,  non-contact,  proximity  sensor 
subsystem  for  monitoring  die  thermally  induced  growth  of  the  spindle  of  a 
machining  center  in  the  X-,  Y-  and  Z-axi$  directions,  plus  the  thermally 
induced  pitch  and  yaw  of  the  spindle, 

(2)  devdopment  of  a  thermal  couple  array  subsystem  to  monitor  and  correlate  the 
thermal  map  of  the  macluning  center  to  the  spindle  growth, 

(3)  development  of  a  personal  computer-based,  user-friendly  data  acquisition  and 
thermal  modeling  system  to  characterize  the  thermal  growth  behavior  of  the 
machining  center  and 

(4)  testing  and  demonstration  of  the  system  in  different  production  environments. 

As  a  result  of  this  program,  an  advanced,  low-cost,  spindle  thermal  growth 
modeling  prototype  ^em  was  developed.  The  system  includes: 

(1)  a  S-ans  capadtance  sensing  and  an  integrated  circuit-based  thermal  sensing 
subsystem  to  monitor  the  spindle  thermal  growth, 

(2)  an  IBM-PC-based,  14-bit,  64-chaimel,  real-time  data  acquisition  subsystem, 
and 

(3)  a  computer-aided  sensor  correlation  (CASC)  software  system  for  sensor 
optimization,  and  modeling. 

The  system  was  fully  tested  and  demonstrated  at  a  number  of  locations  which 
included  the  National  Institute  of  Standards  and  Technology  (NIST),  Pratt  &  Whitney, 
Univeraty  of  North  Carolina,  Olofsson  Corporation,  Ford  Motor  Company  and  The 
Monarch  Machine  Tool  Company. 


2 


■III  I  U.J  j]  iiiMj— 


"1 


Tests  indicated  that  the  ^em  hardware  was  very  reliable  and  very  easy  to  use. 
The  mod^Qg  process  generally  took  less  than  three  days  and  the  modeling  accuracy  was 
about  80%.  This  required  time  was  only  about  one-twentieth  of  the  time  taken  for  most 
modding  efforts,  even  ^en  undertaken  in  a  research  environment. 

These  encouraging  test  results  eventually  led  to  the  development  and  introduction 
of  a  new  product  -  the  THERMAC  (Spindle  Thermal  Modeling  and  Compensation 
System)  ^  API.  Details  of  the  system  and  the  modeling  process  are  described  in  later 
sections  of  this  report. 

2.  PRIOR  WORK 

It  is  generally  agreed  that  the  thermal  expandon  of  a  machine-tool  will 
dramatically  affect  the  machining  accuracy.  After  conducting  an  extensive  study,  Peklenik 
[1]  remarked  that  thermal  errors  could  contribute  as  much  as  70%  of  the  machining  error. 
Therefore,  among  all  the  major  machining  oror  sources,  the  thermally  induced  error  is 
considered  to  be  the  major  problem.  There  are  two  methods  that  are  commonly  used  in 
reducing  the  thermally  induced  errors  on  machine  tools:  error  avoidance  and  error 
compensation. 

Error  avoidance  techniques  are  often  implemented  at  the  machine-design  stage. 
These  indude  designs  to  minimize  the  effects  of  heat  generation  within  the  machine 
structure  and  control  the  gradient  of  and  change  in  environmentally  encountered 
tenq)eratures.  J.B.  Bryan  [2]  of  Lawrence  Livermore  Laboratory  was  successful  in 
achieving  thermal  stability  by  immersing  the  >\diole  machine  in  a  temperature-controlled 
"ofl  bath."  However,  this  method  is  very  expendve  and  impractical  for  production 
maclunes.  A  less  expensive  and  more  practical  method  implemented  by  some  machine- 
tool  builders  on  more  expensive  machines  makes  use  of  a  temperature-controlled  oil 
shower  applied  around  the  spindle  area  where  the  nuyor  thermal  growth  is  typically  found. 
In  generaJ,  this  method  can  help  reduce  the  thermal  growth  error  by  as  much  as  S0%.  Of 
course  the  cost  of  the  machine  is  increased. 

Error  compensation  techniques  are  based  on  a  knowledge  of  the  thermal  growth 
characteristics  of  the  machine  and  use  this  "thermal  model"  in  real-time  to  predict  the 
themud  growth.  A  distal  processor  utilizes  tins  model  to  calculate  the  error  information 
which  is  then  sent  to  the  machine's  CNC  controller  to  compensate  for  the  sensed  error. 
This  error  information  is  often  in  the  form  of  tool  offset  commands. 

By  using  multiple  thermal  sensors  to  measure  the  temperatures  of  the  major  heat 
sources,  and  five  proximity  sensors  to  measure  the  spindle  thermal  growth  (in  the  X-,  Y- 
and  Z-axis  directions,  and  in  pitch  and  yaw),  a  computer-controlled  modeling  system  will 
collect  and  store  the  data  while  the  machine  spindle  is  running.  After  enough  data  is 
coUected,  the  system  will  establish  a  thermal  model  by  analyzing  the  correlation  between 
the  temperatures  and  the  thermal  growth  of  each  axis  of  the  machine.  As  soon  as  the 
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performance  of  the  model  is  verified,  the  model  can  be  used  in  real-time  to  predict  the 
thennal  growths  of  the  machine  by  using  only  the  tempoature  information. 


The  error  compensation  technique  has  raised  a  lot  of  interest  among  machine  tool 
mamifiicturers  and  users.  It  has  the  advantages  of  being  low  cost,  easy  to  implement  and 
can  be  readily  retrofitted  to  existing  machines  without  major  design  modifications.  Real¬ 
time  thermal  growth  compensation  can  be  implemented  using  today's  Mgh  speed  desk  top 
personal  computm  or,  ahematively,  by  using  a  specially  designed  single  board  computer. 

Recent  research  on  thermal  error  compensation  conducted  by  F.  Rudder  and  A. 
Donmez  [3]  of  NIST  involved  the  establishment  of  a  real-time  geometric/thermal  (G-T) 
model  to  characterize  the  thermal  behavior  of  a  vertical  spindle  machining  center  and  a 
turning  center.  The  geometric  error  data  at  various  thermal  states  of  the  machine 
conditions  were  acquired  and  used  to  build  the  G-T  model.  J.S.  Chen  and  J.  Ni  [4]  of 
the  Univerrity  of  hfichigan  (UM)  have  recently  developed  a  time-variant  volumetric  error 
model  which  synthesizes  both  the  geometric  and  thermal  errors  of  a  machining  center. 

In  the  case  of  the  UM  research,  a  21  parameter  kinematic  model  was  used  to 
establish  the  geometric  model  of  the  machine.  An  1 1  parameter  thermal  growth  model 
was  used  to  establish  the  spindle  growth  parallel  to  the  primary  axes  of  the  machine.  All 
together,  a  32  parameter  G-T  model  was  established.  In  both  cases,  tests  were  conducted 
in  a  laboratory  environment  under  dry,  no  load  conditions.  A  laser  interferometer  system 
was  used  for  the  21  geometric  parameter  measurements,  and  a  combination  of  laser  and 
capacitance  measurements  were  used  for  the  1 1  parameter  thermal  growth  measurements. 
The  models  took  2  to  6  months  to  establish,  far  too  long  to  be  practical.  The  accuracy 
improvement  was  about  90%  for  simple  warm-up  and  cool-down  cycles  of  the  machine. 

The  strategy  adopted  by  API  is  somewhat  different  from  that  taken  by  NIST,  UM 
and  others.  Instead  of  establishing  a  G-T  model,  the  API  technique  simplifies  the  entire 
process  by  focusing  on  the  thermal  model  alone.  A  new  model  was  developed  which 
utilizes  d^erential  temperature  data  as  input  parameters.  This  model  has  proven  to  be 
very  effective  in  handling  actual  production  situations  involving  such  things  as  short 
duration  speed  changes,  spindle  stops,  ramps  up  and  down,  etc.,  as  well  as  the  long-term 
warm-up  and  cool-down  cycles.  By  focusing  on  the  thermal  model  alone,  the  entire 
testing  and  modeling  processes  are  drastically  reduced  and  the  accuracy  enhanced. 

To  establish  the  geometric  model,  some  well-known  techniques  and  measuring 
systems  like  the  API  WINNER  2.0  system  can  be  used.  Compared  to  the  independent 
laser  interferometer  measurements  that  were  used  at  NIST  and  UM,  WINNER  utilizes 
alternative  measuring  devices  such  as  a  telescopic  ballbar,  electronic  autocollimator, 
electronic  level  and  others  to  measure  the  geometric  errors  of  the  machine.  In  cases  like 
squareness,  straightness  and  flatness  measurements,  these  alternative  devices  are  much 
better  suited  for  the  job  than  the  laser  interferometer.  They  are  less  costly,  much  easier  to 
setup  and  sometimes  more  versatile  than  a  laser  interferometer  system.  Consequently,  the 
combination  of  simplified  but  advanced  modeling  procedures,  and  a  better  choice  of 
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meiauimg  iiMttumentg  has  atowed  the  API  tediroiue  to  shortened  the  actudn^^ 
process  freun  nxxidu  u>  a  matter  of  a  few  di^. 

3.  TECHNICAL  OBJECTIVE 

The  program  objective  was  to  devdop  an  advanced,  commercial,  easy-to-use, 
inexpensive  j^indle  error  characterization  and  compensation  system  for  production  CNC 
machining  centers. 

To  achieve  this  objective,  it  is  usefel  to  fest  conader  the  spindle  of  a  machine  tool 
to  be  a  rigid  body.  Then  the  spindle's  total  motion  in  a  three-dimensional  machine  work 
zone  can  be  represented  using  6  d^rees-of  -fireedom,  i.e.,  three  translation  motions  along 
the  X,  Y  aiKl  Z  axes  of  the  machine  slides,  plus  pitch,  yaw,  and  roll  motions  around  these 
three  tvthogonal  axes.  Seelngurel.  Since  roll  is  the  turning  axis  of  the  ^indle,  it  is  not  a 
measuring  parameter. 

In  order  to  accomplish  the  program  objective,  the  followings  migor  tasks  were 
pursued; 

(1)  development  of  a  low-cost,  high-performance,  non-contact,  proximity  sensor 
subsystem  for  monitoring  the  X-,  Y-  and  Z-axis  thermal  growth  of  the  spindle 
of  a  machining  center,  phis  the  associated  pitch  and  yaw  growth, 

(2)  devdopment  of  a  thomocouple  array  subsystem  to  monitor  and  correlate  the 
thermal  nuq>  of  the  machining  center  to  the  s^nndle  growth, 

(3)  devdopment  of  a  PC-based,  user-focndly  data  acquisition  and  thermal 
modetog  system  to  characterize  the  thermal  growth  bdiavior  of  the  machining 
center, 

(4)  testing,  demonstration  and  evaluation  of  the  system  on  CNC  madiinii^ 
colters  in  a  production  environment. 

Figure  2  illustrates  the  link  between  the  hardware  and  software  of  the  prototype 

system. 

4.  SENSOR  SYSTEM  DEVELOPMENT 

4.1.  Proximity  Soisor  Selection  &  Development 
Proximity  Sensor  Requiremoits 

Optical  and  ciqiadtance  sensors  are  the  two  most  commonly  used  sensors  for  high- 
precision  non-ccmtact  measuronent  in  industiy.  In  the  early  stage  of  development,  efforts 
were  directed  toward  the  testing  of  various  proxumty  sensors  for  that  ^plication.  Key 
fectors  that  needed  to  be  consideed  were: 

(a)  Measuring  range  -  Typical  ^indle  drift  is  between  ±0.003"  and  ±0.010.” 

Therefore  a  measuring  range  of  at  least  ±0.015"  was  reqiured. 
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(b)  Sensitivity- Sensitivity  ofthe  sensor  should  be  0.000040”  (1/nn)  or  better.  This 
can  cover  a  wide  range  of  machine  accuracy  requirements  and  machine  types,  and 

(c)  Rotaistness  -  Sensor  shcHild  be  fiurtoiy  hardened  and  robust  enough  for  production 
applications. 

Optical  Di^lacement  Sensors: 

The  first  non-contact  sensor  tested  was  an  optical  sensor.  The  sensor  was 
constructed  based  on  a  Hewlett-Packard  HBCS-1 100  optical  device.  Its  measuring 
principle  was  based  on  an  image  focusing  and  defocuang  effect.  See  the  performance 
characteristic  shovm  in  Figure  3.  By  limiting  the  measuring  range  on  either  side  of  the 
focus  point  of  the  sensor,  absolute  di^lacem^t  measurements  could  be  made  with 
reasoiiable  accuracy. 

Based  on  the  tests,  the  following  w^  the  characteristics  of  the  optical  sensor: 
range:  2  mm 

sensitivity:  1  tim  per  millivolt 

voltage  output:  10  volt  dc 

short-term  stability:  better  than  1  millivolt 

The  advantages  of  the  optical  sensor  w&e  that  it  was  easy  to  use,  h  was 
ine^qrensive  and  it  was  based  on  a  compact  design.  The  nuyor  disadvantage  for  this 
^pUcation  was  that  the  senritivity  of  the  sensor  changed  rignificantly  with  different 
surfiu:e  conditions  and  incident  angles.  Conditions  such  as  differait  surfiice  finishes, 
surfiu:e  contamination,  and  surfiice  angle  compared  to  normal  could  alter  the  sensor 
calibration  very  easily. 

Two  optical  sensors  were  initially  tested  on  a  lathe  at  API.  See  Figure  4  for  setup. 
The  sensors  were  mounted  on  an  aluminum  rectangular  fixture  &cing  to  a  measuring  bar. 
Results  indicated  that  while  the  optical  sensors  did  provide  enough  sensitivity  and  range 
for  the  spindle  thermal  error  measurement,  thdr  output  rignals  were  greatly  affected  by 
the  surfiice  conditions,  such  as  contamination,  of  the  measuring  bar.  By  contaminating  the 
bar  surfime  with  a  light  film  of  lubricating  oil,  the  sairiti>dty  of  the  sensor  dropped  by  as 
much  as  40%.  This  indicated  that  optical  sensors  were  undesirable  for  use  in  a  production 
environmoit. 

Capacitance  Prommity  Sensors: 

Ciqiachance  sensors  are  very  precise,  non-contact  devices  but  are  oftoi  very 
delicate.  Most  commerdal  capacitance  sensors  are  not  suitable  for  spindle  thermal  drift 
measurement  because  of : 

(1)  not  having  enough  measuring  range  (typically  below  0.004") 

(2)  having  significant  thermal  drift  characteristics  (most  c^acitance  sensors 
drift  by  as  much  as  40  p  inches  per  "C),  and 

(3)  being  too  costly  for  the  application  (^proximately  $5,000  per  sensor). 
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Tlw  previously  expressed  concerns  about  optical  sensors,  and  the  problons  with 
the  commei^  cq)acitance  sensors  just  mentioned  eventually  led  to  the  development  by 
API  of  a  new,  low-cost,  thermally  stable  and  robust  capacitance  sensor. 

Figure  S  is  the  schematic  of  the  API  capadtance  sensor.  The  basic  differoice 
between  the  commercial  capadtance  gage  and  the  API  sensor  is  that  the  former  uses  a 
bridge  balancing  piindple  while  the  later  uses  a  differential  measurement  prindple.  The 
result  is  that  the  latter  is  much  simpler,  is  less  costly  to  build  and  has  a  greater  operating 
range.  A  drawback  is  that  the  sensitivity  of  the  new  sensor  is  about  ten  times  less  than 
with  the  bridge  versions.  Nevertheless,  the  sensitivity  of  the  new  gage  is  considered 
adequate  for  this  application. 

Output  drift  of  a  capadtance  sensor  is  typically  caused  by  thermal  instability  of  the 
amplifier  drcuits  and  of  the  enwonment.  To  eliminate  the  drift  problem,  a  spedally 
dedgned  compensation  circuit  was  built  into  the  new  sensor.  By  employing  as  a  r^erence 
indde  the  sensor  assembly  a  matching  capacitance  circuit  and  matching  measuring 
element,  the  drift  effect  of  the  primary  sensor  element  is  canceled  out.  See  Figure  6.  As  a 
result,  the  API  capadtance  gauges  experience  only  one-tenth  the  amount  of  drift  that  most 
commercial  capadtance  gauges  exhibit.  V^th  all  ^e  basic  circuits  located  within  the 
gauge  head,  the  API  gauge  can  measure  with  a  much  higher  dynamic  bandwidth  and  with 
better  signal-to-noise  ratio  than  most  commerdal  gauges. 

Figure  7  dqricts  a  drift  test  of  an  API  gauge  using  a  10®C  temperature  variation 
over  a  2  hour  period.  The  gauge  showed  a  thermal  stability  improvement  of  ten  times 
over  that  of  any  of  the  commercial  units. 

Figure  8  shows  the  setup  of  five  of  these  API  capacitance  gauges  when  used  for 
machine  tool  spindle  growth  measurements. 

The  foUo^g  are  the  characteristics  of  the  API  capadtance  gauge: 

Measuring  range:  0.05  mm  to  0.8  mm 

Resolution:  0.10  pm 

Linearity:  better  than  1%  of  full  range 

Thermal  drift:  less  than  100  ppm  per  ®C 

Standoff:  0.05  mm  nominal 

Frequency  response:  20  kHz 

4.2.  Thermal  Sensor  Sdection  &  Development 

After  testing  and  comparing  three  different  types  of  thermal  sensors,  i.e., 
thermocouple,  thermister  and  integrated  circuit  (IC)  thermal  sensors,  the  IC  sensors  were 
selected  for  this  program.  Reasons  for  the  selection  of  IC  sensors  are  that  they  are: 

(a)  more  rugged,  stable  and  accurate  than  thermisters, 
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Figures.  Schonatic  Diagram  of  an  API  Ciq>adtance  Gauge 


Hgure  6.  Block  Diagram  of  a  Ciq;)acitance  Gauge  Circuit 
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(b)  more  suitaMe  for  the  tenqMfsture  mige  and  sensitivity  neecb  anocuted  wi&  ^ 
machme  tool  qiplication  dian  thennocou}^  and 

(c) ea8ytoiitt^ 

The  firflowing  are  the  diaracteristtcs  of  the  IC  senstK  tested: 

Accuracy:  iO.l 

Linearity:  0.1  **C 

Operating  range:  -25  **€  to  +105 

Suppfy  voltage:  +4  Volt  to  +30  Vdt  dc 

A  unique  magnetic  base  design  was  developed  to  house  each  thermal  sensor. 

Figure  9  shows  an  IC  thetmal  sensor  eorf)edded  inside  a  magnet.  To  use  die  sensor,  the 
opeatxx  simpfy  places  the  sensor  against  a  fisrous  madune  tool  surfiice  and  the  magnetic 
^ce  autcnnaticdty  Ixrfds  the  sensor  in  place. 

43.  Sensor  InterfiKe  Controler  Development 

A  high-perfixmance  computer  inter&ce  controller  was  developed  for  the  sensor 
systenL  It  was  incorporates  a  14^  D/A  ccmverter,  paralld  intetfiioe,  and  an  internipt- 
driven  DMA  (direct  memory  access)  interfiice  driver.  Using  a  multiplexing  tedmique,  the 
contixrfler  is  c^iable  of  handling  as  mai^  as  27  thermal  sensors  and  5  capacitance  sensors 
siinultaneously.  Four  control  lines  fiom  the  conqiuter  are  used  for  the  dumnd  sdect. 
Figure  lOisasdiematicdiagramofthesyrtemstowiqgthecomnrfler. 

4.4.  User-FMendly  Software  Dcvdopmenf 

User-fiiendly  interacdve  software  was  devdoped  for  the  sensor  interfiice,  system 
setup,  data  acquisition  and  analysis  fiinctions.  Figure  1 1  illustrates  the  software  structure. 
There  are  fixir  interactive  diqdi^  screens  to  fiidlitate  use  of  the  system.  Thesescreens 
fiidlitate  sensor  setup,  data  coU^on,  programming  and  data  an^^sb.  An  exanqrfe  of  the 
sensor  setup  screen  is  presented  in  Rgure  12.  Senstxs  1, 2, 3, 4  a^  5  represent  the  XI, 
Yl,Z,X2and^ouqxitsofthefivecapadtancegau^.  They  monitm’ the  distance  to 
thesmfooeof  the  test  mandrel  i^di  is  a  polished  1”  diameter  x  8"  long  sted  rod 
mounted  in  the  ^rfndle  of  the  madiine.  During  setup,  the  readouts  from  all  five  gauges  are 
continuously  di^lsyed  cm  the  screen,  thus  9ving  visual  fisedbadt  to  the  user  as  he/she 
acgusts  the  sensors  into  position.  The  total  required  setup  tune  is  about  5  rniinites. 

The  software  program  {uovides  the  firflowing  user  sdectalrfe  fimctkxis  on  the  main 

soeen: 

Data  Collection 

One  can  sdect  the  samplii^  interval,  test  runmng  time  and  the  output  filename. 
During  data  ctrflection,  the  diqili^  will  show  the  ouqxits  of  all  five  oqMdtance 
gages  and  the  user-sdected  thern^  sensors.  To  diininate  the  infiuence  of  the 
bundle  vibration  caused  by  off-coiter  daiddng  of  the  mandrd,  the  ouqxits  of  the 
capacitance  gauges  are  read  and  averted  500  times  over  10  spindle  revolutions 
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Electronics  Controller 
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Figure  10.  Schematic  Diagram  of  Electronics  Controller 


(iurii^  each  sunpliiig  interval.  The  system  is  capaUe  of  detecting  the  actual  rpm 
of  the  madiine  q;undle  and  sets  the  sampiii^  fie^ency  according. 


AnalyMK 

At  the  conclusion  of  the  data  collection  function,  the  results  will  be  di^layed  and 
analyzed  uang  the  analym  screen.  The  data  points  are  computed  to  yield  the  X, 

Y,  Z,  (Mtch  (using  a  conibination  of  the  XI  and  X2  sensor  ou^Hits)  and  yaw(  using 
a  ccMninnation  of  the  Y1  and  Y2  sensor  ou^nits)  movements  of  the  ^indle.  Hard 
001^  of  the  results  can  be  obtained  throuj^  use  of  dot  matrix  or  huw  printers. 

The  ^em  was  tested  at  many  pla'~/es,  and  comments  were  collected.  Overall,  tb» 
hardware  and  scrftware  received  fiivorable  comments  from  its  users.  Suggestions  for 
in^Kovements  were  noted  and  these  ideas  were  incorporated  in  subsequent  vernons  of  the 
Spindle  Thermal  Characterization  and  Compensation  System. 

5.  THERMAL  MODELING  AND  THE  COMPUTER  AIDED  SENSOR 
CORRELATION  SYSTEM  (CACS) 

5.1.  Analytical  vs.  Parametric  Thermal  Modeling  Approaches 

Two  dasacal  approaches  to  modeling  the  thermal  distortion  of  machining  centers 
are  tiie  analyticai  (tinhe  element)  approach,  and  the  parametric  (experimortal)  approach. 

In  the  analytical  approach,  it  is  assumed  that  the  machine  structure  is  homogeneous  and 
that  the  h^  disapation  throughout  the  structure  follows  the  laws  of  thermal  conductivity. 
These  assumptions  are  then  used  to  establi^  an  ideal  analytical  model  (i.e.,  finite  demem 
modd)  of  the  machine.  Users  of  this  model  first  identify  foe  locations  of  foe  heat  sources 
and  estimate  the  amount  of  heat  dissipated.  Then  foe  model  predicts  foe  total  thermal 
distortion  [5]. 

Unfortunately,  researchers  [6]  have  found  that  foe  boundary  conditions  of  two 
mating  surfoces,  which  often  greatly  affect  the  predicted  thermal  transmissivity,  are  very 
difficult  to  estimate  with  good  accuracy.  In  addition,  it  is  also  difficult  to  estimate  the 
amoum  of  heat  generated  and  dissipat^  without  making  actual  measurements  on  foe 
madune.  The  lade  of  predse  knowledge  of  the  rrumy  boundary  conditions  and  foe 
difficulties  assodated  with  heat  disdpation  foctors  often  cause  skeptidsm  about  foe  results 
obtained  from  this  analytical  approach.  Consequently,  foe  analytical  approach  is  used 
iiutnfy  for  trend  predictions  nfoer  than  for  qualitative  predictions. 

The  work  presented  in  this  report  has,  foerdbre,  been  based  on  foe  parametric 
modding  approach  which  requires  no  knowl^ge  of  foe  machine  configuration,  structure 
and  boundary  conditions.  The  approach  uses  actual  measurements  of  foe  thmmal  and 
growth  effects  of  foe  machine  spmdle.  Throu^  a  sequence  of  well  designed  tests  and 
ejqretiments,  the  approach  correlates  foe  thermal  and  spindle  growth  information  and 
establishes  a  mathematical  (or  parametric)  model  of  the  machine.  Once  verified,  the  model 


16 


can  dnen  be  used  in  real-time  unth  thermal  meanirements  to  predict  the  thermal  growth  of 
the  machine  and  movement  of  the  spindle. 

The  machine  must  be  exposed  to  a  variety  of  operating  conditions  and 
environmental  conditions  during  the  modeling  process  in  order  to  build  up  a  robust  model. 
In  addition,  the  accuracy  of  the  model  is  very  dependent  upon  the  performance  of  the 
measuring  instruments.  Fortunately,  all  these  are  easily  controlled  with  greater  certainty 
than  would  be  possible  with  the  analytical  approach.  Because  of  this,  the  parametric 
tq)proadi  has  become  more  widely  accepted  by  industry.  Based  on  the  results  of  this 
work,  the  accuracy  achieved  by  the  parametric  approach  was  well  above  80%,  several 
times  more  accurate  than  could  be  expected  from  the  aiudytical  approach. 

5.2.  Spindle  Thennal  Growth  Measurement 

Extensive  tests  of  the  spindle  thermal  growth  were  carried  out  on  various  types  of 
marhitiing  centers  during  and  t^er  the  development  of  the  initial  sensor  system.  The  tested 
machining  centers  included: 

•  an  Autonumeric  MVC-10  vertical  spindle  CNC  machining  center, 

•  a  Monarch  VMC-75  vertical  spindle  CNC  machining  center, 

•  a  Sundstrand  series  20  Omnimill  horizontal  spindle  CNC  machining  center, 

•  a  Producto  A-1738  vertical  spindle  CNC  machining  center,  and 

•  an  Olofrson  dual-vertical  spindle  turning  center. 

Except  for  the  Sundstrand  machine,  these  machining  centers  are  classified  as  "C 
type  vertical  machining  centers.  It  was  discovered  that  the  major  direction  of  spindle 
growth  of  the  "C"  type  machines  is  along  the  Z-axis.  This  growth  is  about  three  to  four 
times  the  growth  in  the  direction  of  the  X-  and  Y-axes.  As  a  consequence,  most  of  the 
data  shown  here  are  for  growth  along  the  Z-axis. 

The  thennal  tests  for  each  machining  center  typically  consisted  of: 

•  Warm-up  run  -  Running  the  spindle  afrer  the  spindle  was  shut  down  overnight, 

•  Run  with  various  spindle  speeds  -  Running  the  spindle  at  different  speeds  after 
the  nuiclune  is  fully  warmed-up, 

•  Production  cycle  run  -  Simulating  the  production  cycle,  including  starting  and 
stopping  the  spindle  frequently, 

•  Cool-down  run  -  Monitoring  the  spindle  movement  after  the  spindle  was 
stopped. 

Figure  13  shows  two  observations  of  the  Z-axis  warm-up  spindle  growth  of  the 
Autonumeric  machine.  Spindle  speeds  were  1000  and  2000  rpm.  Figure  14  shows  a 
typical  warm-up  and  cool-down  test  of  the  Autonumeric  machine  in  both  Y-  and  Z-axes. 

It  is  common  to  find  that  the  warm-up  growth  is  more  rapid  than  the  cool-down 
growth/contraction  since  the  former  is  a  result  of  "forced"  heat  transfer  while  the  latter  is  a 
result  of  "natural"  heat  conduction  and  convection. 
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Figure  13.  Effect  of  Spindle  Speed  Variation 


EFFECT  OF  TOMPERATURE  HBTORT 


Figure  14.  Effect  of  Temperature  Ifistoiy 
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Environiiiemal  eflSects  can  also  significantly  affect  tnaciune  growth.  This  is  shown 
in  Hgure  IS.  The  Sundstnmd  machining  center  was  uikler  a  routine  warm-up  test  when 
suddenly,  at  about  180  minutes  into  the  test,  a  large  bay  door  about  200  feet  away  from 
the  machine  was  op^ied.  The  door  opening  showed  chUled  air  from  outside  to  rush  in 
and  upset  the  mad^  growth  pattern.  As  can  be  seen  from  the  figure,  most  tonpaature 
gauges  re^nded  to  the  suddoi  drop  in  tenqierature  at  180  minutes.  In  the  meantime,  the 
capadtance  sensors  also  picked  up  rapid  changes  in  the  spindle  portion  (large  dips  shown 
by  SI,  S2,  S3  and  S4).  Shortly  after  the  door  was  closed,  both  the  temperature  and  the 
^nndle  portions  returned  to  thdr  previous  status.  Although  there  was  insufficient  data  to 
esqplain  why  the  Z-axis  (monitored  by  SS)  didn't  seem  to  change  as  much  as  the  other  axes, 
the  effect  of  the  environmental  tonperature  fluctuation  on  the  machine  growth  was  vividly 
demonstrated  by  this  incident. 

53  Thermal  Modeling 

It  was  discovered  during  the  many  tests  on  different  machines  that  the  heat  source 
that  seoned  to  influence  spindle  growth  ffie  most  was  the  spindle  bearings.  Most  modem 
machines  have  thdr  spindle  motor  isolated  from  the  main  structure  of  the  machine  and, 
therefore,  the  heat  generated  by  the  motor  has  less  effect  than  the  heat  generated  by  the 
^indle  bearings. 

Our  findings  also  indicated  that  enviroim^tal  effects,  residual  thermal  effects  in 
the  spindle  and  column,  such  as  embedded  heat  accumulated  during  the  previous 
operations  gradually  dissipating  from  the  ma^,  and  current  operating  conditions  can  also 
greatly  affect  the  spindle  growth. 

Due  to  the  difference  in  mass  between  the  spindle  and  column,  the  spindle 
normally  heats  up  and  cools  down  much  quicker  thm  the  column.  This  makes  thermal 
modeling  voy  complicated,  especially  when  the  machining  center  is  used  m  a  production 
environment.  In  a  typical  production  environment,  the  machine  is  often  subjected  to  the 
influence  of  more  conqrlex  thermal  growth  processes,  other  than  simple  warm-up  and 
cool-down  cycles,  and  sometimes  uncontrolled  environmootal  conditions.  These  are  often 
the  causes  of  fiiilure  of  the  conventional  parametric  modeling  methods. 

Based  on  these  findings,  API  believed  that  the  modeling  had  to  be  performed  in  an 
environment  that  the  machine  was  intended  for.  As  a  further  consideration,  the  residual 
thomal  effect  and  the  speed  sensitive  effect  which  could  be  measured  using  differential 
thomal  measurements  along  the  machine  spindle  and  column  should  be  an  included 
parameter  in  any  newly  developed  thermal  modeling  techmque. 

Subsequently,  a  new  multiple  regresaon,  differential  thermal  parameter  modd  was 
developed. 
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Hie  general  fbnn  of  the  model  for  a  single  axis  is: 

Etm  =  Aom  +  AlmATit  +  AlmAJ^  +  AlmhTid  -f  A*mhJ^ 

Where: 

£ta  s  ihermcd  growth  errors  predicted^  where  m  is  mtchine  axis,  i.e.  X,Y,Z... 
Aom.  Aim.  Aim.  Aim.  AAm:  coefficieTits  of  model, 

ATit  —  til  — tit 

tit  is  the  environment  temperature 
tttis  the  selected  thermc  J  setter, 

ATu  =  (/«  —  /# ) 

IJ  =  thermal  sensor  number, 
i  =  time  stamp. 

The  equation  calls  for  the  use  of  multiple  thermal  measurements  for  all  key  heat 
sources  of  the  machine.  In  order  to  ensure  that  all  key  locations  are  measured,  it  is 
common  to  find  a  large  number  of  redundant  thermal  sensors  used  initially.  At  NIST  [3], 
SO  thermal  sensors  were  initially  applied  to  the  machine.  After  days  or  sometimes  weeks 
of  thermal  testing,  a  special  computer  algorithm  was  used  to  correlate  the  thermal  and 
displacement  sensor  outputs.  The  results  were  then  used  to  identify  the  sensors  that  had  a 
hi^  correlation  with  machine  growth,  and  to  eliminate  those  that  luul  a  low  correlation. 
The  elimination  decision  process  was  veiy  lengthy  and  partly  done  manually.  The  final 
number  of  sensors  select^  usually  ranged  firom  four  to  five. 

At  the  Univertity  of  hfichigan  [5],  approximately  50  thermal  sensors  were  used  in 
the  initial  stages  of  the  test.  However,  the  UM  sensor  elimination  process  was  done  in 
two  steps.  The  first  involved  the  use  of  a  similar  correlation  program  to  identify  the  key 
sensors.  The  second  step  was  a  stepwise  regression  analysis  which  identified  the  cross- 
correlation  results  between  thermal  sensors  before  they  were  further  reduced.  Both 
procedures  were  tedious  and  the  results  often  relied  on  the  judgment  of  the  operator. 

To  be  able  to  automatically  and  effectively  reduce  the  number  of  the  thermal 
sensors  on  the  machine  is  a  necessary  part  of  any  commercially  viable  spindle-thermal 
error  analytis  and  compensation  system.  Users  of  the  system  should  not  find  the  modeling 
process  too  cumbersome  to  perform  nor  the  analysis  too  theoretical  to  understand. 
Consequently,  as  a  part  of  the  analysis  and  compensation  system,  API  developed  a  user- 
fiiendly  computer-aided  sensor  correlation  (CASC)  algorithm  to  assist  in  the  optimization 
(minimization)  of  the  number  of  thermal  sensors.  CASC  is  a  real-time  user-interface 
program  designed  to  help  the  user  to  identify  the  key  locations  of  the  thermal  sources  on 
the  machine  as  wdl  as  optimize  the  number  of  thermal  sensors  required.  The  initial 
number  of  thermal  sensors  used  in  the  API  system  was  20  and  the  final  number,  based  on 
our  maiQ^  tests,  was  only  4  or  S. 
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5.4  Principle  of  Sensor  Corrdation 

Correlation  analysis  is  a  statistical  method  used  to  reflect  the  degree  of  the  linear 
relationship  between  two  variables.  The  correlation  coefficient  denoted  by  y  is  defined  as: 


^(xi-xxyi-  y) 

Where: 

Xi^i  =  values  of  variables  x  and  y, 

x,y^  mean  value  of  Xiondyi, 

Y  =  correlation  coefficient  between  x  and  y. 

The  value  of  the  correlation  coefficient  y  is  zero  if  x  and  y  are  random  occurrences 
and  are  unrelated  to  each  other.  The  value  of  y  is  1  if  they  are  totally  related,  i.e.,  if  y  = 
f(x).  See  Figure  16. 

In  the  case  of  C  ASC,  let  y  be  the  correlation  coefficient  between  the  thermal 
sensors  and  the  spindle-axis  growths.  Then  y  has  the  form: 

Where: 

=  values  of  thermal  sensor  ti,  I  is  sensor  number  from  1  to  20, 
ti  =  mean  value  of  tu, 

dk  =  values  of  ccpacitance  sensor  dk,  k  is  sensor  number  from  1  to  5, 

dk  =  mean  value  of  d*, 

ytk  =  correlation  coefficient  between  ti  and  dk, 

i  =  time  stamp. 

CASC  calculates  the  correlation  coeffidents  between  the  thermal  drift  as  measured 
by  each  capacitance  gauge  and  the  temperature  as  measured  by  each  thermal  sensor.  It 
then  prioritizes  the  thermal  sensors  (from  high  to  low)  based  on  the  y  values.  Figure  1? 
iUustrates  the  results  of  the  Z-axis  correlation,  temperature  measurements  against  Z-axis 
growth,  produced  by  CASC  on  the  Monarch  machining  center.  The  correlation  values  are 
shown  as  bars  above  the  horizontal  line  while  the  bottom  half  of  the  graph  shows  the 
actual  temperature  range  measurements.  The  sensors  were  prioritized  based  on  thdr  y 
value,  and  plotted  with  the  highest  value  sensor  on  the  tight.  As  can  be  seen,  thermal 
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Figure  17  CASC  System  Graphics  Display 


sensor  #13,  ^ch  had  the  highest  tenqierature  increase,  with  +1S  X,  also  had  the  hi^iest 
Y  value.  It  was  located  at  the  spindle  nose,  closest  to  the  lower  spindle  bearing.  CASC 
showed  that  sensor  #13  had  the  highest  correlation  to  the  Z-axis  growth  of  the  noachine 
spindle.  Rgure  18  illustrates  the  various  thennal  sensor  locations  on  the  machine. 

It  is  important  to  realize  that  higher  temperature  measurements  do  not 
automatically  imply  higher  y  values.  As  is  illustrated  back  in  Figure  17,  sensor  #7 
experienced  less  tempmture  increase,  but  shown  to  have  a  higher  y  value  than  that  of 
#1.  This  implies  that  the  CASC  provides  a  measure  of  quality,  not  quantity. 

Automatic  sensor  sdection  using  purdy  y  values  could  be  misleading  since  the 
locations  of  the  sensors  on  the  machine  would  not  be  a  &ctor  in  the  evaluation  process.  In 
other  words,  there  could  be  sensors  that  were  located  at  places  where  they  were  simply 
measuring  the  same  heat  source,  and  both  would  present  a  high  y  value.  Such  a  pair  of 
sensors  would  be  #13  and  #10.  To  prevent  this  condition,  it  is  necessary  for  CASC  to 
establish  a  second  correlation  process  to  ch^k  the  correlation  between  thermal  sensors. 
This  process  is  referred  to  as  autocorrelation,  in  other  words,  corrdation  between  the 
same  type  of  variable,  in  this  case  thermal  variation. 

For  this  purpose  the  autocorrelation  coeffident  a  is  defined  as; 

Vs  (*'"-* 


Where; 

ttt  =  values  of  thermal  sensor  ti,  I  is  sensor  number  from  1  to  20, 
ti  =  mean  value  of  tu, 

tik  =  values  of  thermal  sensor  tk,  k  is  sensor  number  I  <k^  20, 
ft  =  mean  value  of  tk, 

cat  =  autocorrelation  coefficient  between  ti  and  tk, 
i  =  time  stamp. 


If  two  thermal  sensors  monitor  roughly  the  same  heat  source,  a  should  be  close  to 
unity.  Based  on  &q)erience.  we  selected  0.97  as  the  threshold  value  for  a.  In  other 
words,  if  two  or  more  sensors  had  an  a  equal  to  or  larger  than  0.97,  then  only  one  would 
besdected.  The  remaining  ones  would  be  eliminated  in  the  modeling.  As  an  example,  the 
autocorrdation  coeffident  of  thermal  sensors  #13  and  #10  was  0.98S.  As  a  consequence, 
we  would  eliminate  dther  #10  or  #13  fi-om  our  final  setup. 

In  sununary,  the  operating  procedures  of  CASC  are; 

(a)  prioritize  the  thermal  sensors  based  on  the  thermal-to-displacement 
correlation  coefficient  value,  y. 
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(c)  eHminate  the  redundant  thennal  tensors  that  have  a  vahies  equal  to  or  greater 
than  0.97. 

Conceivably  the  number  of  thennal  seosma  should  be  increased  in  order  to 
provide  a  robust  enou|^  model.  This  decision  can  be  made  by  looking  at  the  residual 
value  Rofther^ression  analysis.  Note  that  the  value  of  R  is  unity  when  a  perfect  modd 
is  established.  An  iterative  modeling  process  can  be  used  to  ensure  that  the  modd  is 
optimized.  This  is  evidenced  by  having  a  value  for  R  that  is  close  to  unity.  Typicallya 
seccHid  or  sifosequent  modd,  using  the  one  or  two  additional  sensors  lun^  the  next 
higgler  Y  values,  are  generated.  Then  the  R  values  ofthenoodds  are  compared.  Ifthe 
value  ofRinqnoves  by  more  than  1%  over  the  ixevious  noodd,  the  iterative  process  will 
continue.  It  is  our  experience  that  the  qitimizing  process  seldom  requires  more  than  five 
iteratirms.  As  a  exanfie  of  how  the  process  worlc^  the  initial  CASC  modd  may  suggest 
the  use  of  only  two  thermal  sensors,  such  as  one  for  the  table  and  one  for  the  ^indle. 

After  the  iterative  process,  the  optimal  modd  may  require  four  sensors,  namdy  one  for  the 
taUe,  one  ftff  the  spindle  and  one  differential  pur.  Since  the  process  is  conqdetdy 
automatic  and  onfy  deals  with  data  already  stored  in  the  system,  the  whole  {Mocedure 
usually  doesn't  take  more  than  five  minutes. 

6.  TESTING  AND  EVALUATION  OF  THE  THERMAL  MODELING  SYSTEM 

Most  of  the  system  testify  and  evaluation  was  conducted  at  the  program 
coUaborator's  she,  at  The  Monarch  N^Kfaine  Tool  Compare  fiidlity  in  Cortland,  New 
York.  Some  of  the  testing  and  evaluation  was  conducted  at  the  Oloftson  Corporation 
plant  in  Lansing,  Nfichigan. 

Figure  18,  mentioned  previously,  shows  the  sensor  configuration  on  the  Monarch 
VMC-4SBmachming  center  used  fix  the  testing  and  evaluation.  The  machine  was 
equi|^)ed  with  a  qrindle  temperature  ccmtroller,  or  chiller,  capable  of  ddivering 
temperature  controlled  oA  directly  to  the  q;>indle  bearings  and  to  the  axis  slides  at  a 
tenqMrature  hdd  to  within  ±  2  *X).  The  maximum  q)eed  ofthe  yindle  was  10,000  ipm. 

One  would  think  that  the  indusion  of  the  duller  in  the  machine  setup  might  have 
c(Hiq)ticatod  the  moiling  and  compensation  process.  We  fixmd,  however,  that  vdien  the 
^an^  speed  was  above  3000  rpm  the  duller  became  very  ineffective.  Also,  in  order  to 
demcHistrate  the  full  effect  of  the  modding  syston,  most  of  the  tests  were  focused  on 
operations  at  speeds  of 3000  rpm  or  higher. 

The  whole  process  ofmodding  and  evaluation  was  conducted  in  three  stqis.  The 
first  stq>  was  data  collection  <»  a  vertical  machining  center  at  The  Monarch  Maclune  Tool 
Conqpany.  The  seermd  was  data  processing  and  modding  at  API.  And  the  third  was 
modd  evaluation  on  the  same  machinii^  center  back  at  Monarch. 
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With  the  system  initially  tested  at  API  and  several  production  fiidlities,  it  was  then 
taken  to  Monarch  for  the  first  step  of  data  collection  on  the  vertical  machining  cento’.  For 
3  days,  large  amounts  of  data  was  collected  which  included  periods  of  warm-up,  cool¬ 
down,  random  ^>eed  tests,  short-duration  high-speed  tests,  and  short-duration  zero  speed 
tests.  This  was  meant  to  simulate  all  the  different  ^)eeds  and  operating  conditions  that 
might  be  e3q)oienced  on  the  machining  center.  The  spindle  was  found  to  have  four  to  five 
times  more  growth  in  the  Z-axis  when  compared  to  the  growth  along  the  X-  and  Y-axes. 
Thus,  for  illustrative  purposes,  the  modeling  process  was  focused  on  the  Z-axis  growth. 

Next,  the  data  were  brought  back  to  API  and  carefully  analyzed.  The  process  of 
sensor  elimination  and  model  optimization  was  performed  and  results  were  compared  to 
what  would  be  obtained  firom  the  conventional  modeling  technique,  in  other  words, 
without  using  different  sensor  information.  It  was  discovered  that  the  new  model  had  an 
average  improvement  of  approximately  15%  over  the  conventional  technique.  In 
particular,  the  new  model  was  able  to  predict  the  growth  more  precisely  during  random 
speed  operations.  See  Figure  19.  In  this  particular  case,  CASC  had  identified  thermal 
sensors  13, 19,  7  and  4,  and  differential  pairs  13/17  and  10/13  as  being  appropriate  for 
incluaon  in  the  Z-axis  model. 

After  the  analysis  and  modeling  step  was  completed,  the  system  was  brought  back 
to  Monarch  and  set  up  on  the  same  machine  for  the  third  step  of  the  work.  The  machine 
was  once  again  programmed  to  perform  ail  the  same  warm-up,  cool-down  and  random 
speed  cycles  while  the  system  again  collected  data.  However,  this  time,  instead  of  using 
the  data  to  establish  a  model,  they  were  used  to  make  a  comparison  with  the  results 
obtained  fi-om  the  model  previously  established.  Figures  20  through  24  show  the  results 
of  the  comparison  at  spe^  of 3,000, 6,000, 7,000  and  9,000  rpm.  As  can  be  seen,  the 
thermal  growth  was  about  60  pm  with  the  nuudmum  error  caused  by  the  model  being 
7pm.  This  was  an  improvement  of  88%! 

Figures  25  and  26  show  results  of  the  comparison  under  conditions  of  random 
speed  operation.  In  both  situations  shown,  the  machine  was  fully  "warmed-up"  using 
continuous  operation  for  more  than  8  hours.  However,  as  can  be  seen,  even  a  so  called 
"warmed-up”  machine  can  experience  large  thermal  growth  fluctuations  simply  due  to 
varying  operating  speeds.  Figure  25  emulated  a  typical  "on-machine  gauging”  cycle,  a 
modem  manu&cturing  process  that  has  generated  tremendous  interest  in  the 
manu&chiring  community.  The  steep  downward  trends  on  the  gr^h  were  created  when 
the  spindle  was  suddenly  run  at  high  speed.  The  steep  upward  trends  were  created  when 
the  qnndle  was  suddenly  stopped,  a  condition  that  would  simulate  the  start  of  an  on- 
machine  part  gauging  operation.  This  cleariy  illustrates  the  foUacy  associated  with  many 
on-machine  gauging  operations  being  performed  day-in  and  day-out  in  many 
manu&cturing  focilities,  namely,  that  such  measurements  consistently  provide  good  data. 
To  the  contrary,  the  machine  geometry  can  significantly  change  during  the  course  of  an 
on-machine  gauging  operation!  Figure  26  also  illustrates  the  limitation  or  uncertainty  of 
the  presoit  spindle  thermal  modeling  system.  The  best  achievable  accuracy  is  about  ±2.5 
pm  in  a  random  speed  operation. 
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Figure  21.  Thennal  Groivdi  Prediction  in  Z-axis  (60(X)  rpm) 
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Figure  22.  Thermal  Growth  Prediction  in  Z-axis  (6(X)0  rpm) 
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I4gure23.  Thermal  Growth  Prediction  in  Z-axis  (7000  rpm) 
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Kgure24.  Thermal  Growth  Prediction  in  Z-axis  (9000  rpm) 
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Ingure2S.  Thermal  Growth  M(xieling  in  Z-axis 
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Figure  26.  Thermal  Growth  Prediction  in  Z-axis  (random  speed) 


Thermal  growth  tests  of  the  spindle  at  different  Z-axis  locations  were  also 
performed.  It  was  generally  expected  that  at  different  spindle  positions,  the  thermal 
growth  characteristics  might  be  different.  To  verify  this,  the  spindle  was  raised  by  14 
inches,  almost  the  foil  Z-axis  column  travel,  and  similar  spindle  test  runs  were  performed. 
However,  as  can  be  seen  from  Figure  27,  the  difference  in  thermal  growth  characteristics 
were  not  significant  at  all.  This  suggests  that  a  uniform  model  can  be  successfully  applied 
for  this  Monarch  machine.  If  significant  differences  were  observed,  individual  models 
would  have  to  be  established  at  different  positions  of  the  machine.  However,  based  on  our 
testing,  this  would  only  apply  on  some  very  large  machining  cento^s  and  the  improvement 
may  not  be  significant  enough  to  warrant  the  lengthy  procedure. 

Another  foctor  that  might  have  influence  on  the  thermal  growth  behavior  of  a 
machining  center  is  the  heat  generated  by  the  axis  drive  systems,  especially  when  the 
carriages  are  subjected  to  near  continuous  rapid  traverse  rates.  Tests  to  determine  such 
effects  were  carried  out  on  the  Monarch  machine  and  the  Olofsson  machine,  however,  no 
significant  changes  to  the  thermal  patterns  were  observed.  See  Figure  28. 

7.  CONCLUSIONS  AND  DISCUSSIONS 

The  following  statements  summarize  the  major  achievements  of  this  work: 

•  A  low-cost,  prototype  spindle  thermal  characterization  and  compensation  system  for 
CNC  machining  centers  has  been  successfully  developed.  The  system  consists  of  five 
high-performance,  non-contact  capacitance  sensors  mounted  in  an  appropriate  fixture, 
a  mandrel  for  mounting  in  the  spindle  of  the  machine,  a  20-chaimei  thermal  measuring 
unit,  a  series  of  specially  designed  thermal  sensors  capable  of  being  magnetically 
mounted  on  the  machine,  a  PC-based  digital  processor  and  sensor  interface,  and  a 
suite  of  user-fiiendly  thermal  modeling  and  sensor  optimization  software  programs, 
collectively  referred  to  as  CASC. 

•  The  prototype  system  has  been  successfully  tested  and  implemented  on  machining 
centers  in  the  plants  of  several  major  U.S.  manufacturers.  The  projected  accuracy  of 
improvement  is  better  than  80%. 

•  The  entire  modeling  process  usually  takes  less  than  3  days.  This  is  significantly  less 
time  than  that  taken  by  any  other  modeling  technique.  As  a  consequence,  the  use  of 
this  system  on  a  machine  in  a  production  environment  is  very  practical. 

•  The  hardware  and  software  developed  during  tins  program  have  been  successfully 
commercialized  by  API.  Currently  API  is  marketing  this  system  under  the  trade  name 
THERMAC.  THERMAC  is  the  first  low  cost,  high  performance  spindle  thermal 
modeling  and  compensation  system  designed  specifically  for  use  on  machines  in  a 
production  environment.  Of  course  it  can  be  used  on  machines  in  a  laboratory  or  tool 
room  setting  as  well. 


Figure  27.  Thermal  Growth  in  Z-axis  at  Two  Different  Locations 
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Figure  28.  Thermal  Test  Result  with  Slide  Movement 


As  noted,  the  hardware  and  software  of  the  prototype  spindle  thermal  modeling 
system  was  fully  tested  and  implemented  on  several  different  machining  centers.  Most  of 
the  tests  were  performed  in  production  environments.  Some  were  done  in  laboratory 
environments.  Each  test  typically  ran  continuously  for  3  to  S  days.  The  hardware  proved 
to  be  very  robust,  rdiable,  accurate,  easy  to  setup  and  the  software  was  shown  to  be  very 
user-friendly.  Besides  the  tests  done  at  Monarch  and  Olofsson,  some  of  the  initial  tests 
were  also  carried  out  by  engineers  from  Pratt  &  Whitney  and  Ford  Motor  Company. 
Thanks  is  offered  to  all  those  involved  at  each  of  these  facilities  for  their  valuable  input 
and  feedback. 

Based  on  our  testing,  and  for  a  given  machine  structure  and  design,  there  are  three 
sets  of  Actors  that  dominate  the  thermal  growth  pattern  of  a  machining  center.  They  are: 
(a)  the  envirotunental  factors,  (b)  the  residual  heat  factors,  and  (c)  the  operating  condition 
factors  such  as  speed,  operating  cycle,  etc.  As  was  vividly  illustrated  in  Figure  15,  ihe 
machine  geometry  can  be  significantly  distorted  due  to  any  change  in  the  ambient  or 
environmental  temperature.  Additionally,  residual  heat  inside  the  machine  column  and 
spindle  can  cause  non-linear  growth  of  the  spindle,  and  thus  make  the  prediction  and 
modeling  of  the  thermal  characteristics  of  the  machine  very  complicated.  The  machine's 
current  and  previous  operating  conditions  also  significantly  affect  the  thermal  growth 
pattern.  Different  spindle  speeds  cause  different  rates  of  growth.  Different  initial  thermal 
conditions  also  affect  the  growth  patterns,  even  though  the  machine  is  subjected  to  the 
vfty  same  operating  cycles. 

The  numerous  tests  of  spindle  thermal  growth  that  were  conducted  and  the 
numerous  in-depth  analyses  that  followed,  led  to  the  introduction  under  this  program  of  an 
expanded  thermal  modeling  technique  that  includes  the  use  of  differential  thermal  sensors. 
The  use  of  differential  sensors  was  found  useful  in  overcoming  uncertainties  caused  by  the 
three  sets  of  frictors  noted  in  the  preceding  paragraph.  Results  showed  that  the  new  model 
was  able  to  improve  the  modeling  accuracy  by  15%,  even  when  applied  in  a  production 
en\dronment. 

Actual  machine  tool  compensation  for  thermal  growth,  using  the  modeling 
technique  developed  here,  was  not  implemented  and  was  not  a  part  of  the  scope  of  this 
project.  However,  the  accuracy  of  this  modeling  technique  was  verified  by  testing  the 
predicted  results  against  actual  measurements.  Tests  showed  that  the  accuracy 
improvement  of  tlus  spindle  thermal  modeling  system  is  better  than  80%.  Because  of  the 
uncertainty  of  the  thermal  growth  of  machine  structures,  the  accuracy  limit  of  this 
technique  is  about  5  |xm.  In  the  later  stages  of  this  work,  particularly  relating  to  the 
measurements  taken  at  Ford  and  at  Pratt  &  Whitney,  the  entire  modeling  process  took 
only  3  days,  including  the  taking  of  the  geometric  measurements.  Real-time  thermal  error 
compensation  can  be  implemented  by  interfacing  the  prototype  control  unit  through  a 
communications  link  to  the  CNC  controller.  By  sending  the  computed  compensation 
value  to  the  CNC  controller,  the  controller  can  offset  the  tool  position  and/or  orientation 
to  compensation  for  spindle  growth. 
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